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G\ [ ABSTRACT 

We present high-precision photometry of two transit events of the extrasolar planetary 
system WASP-5, obtained with the Danish 1.54 m telescope at ESO La Silla. In order 
to minimise both random and flat-helding errors, we defocussed the telescope so its 
CO ' point spread function approximated an annulus of diameter 40 pixels (16"). Data 

I reduction was undertaken using standard aperture photometry plus an algorithm for 

0^ . optimally combining the ensemble of comparison stars. The resulting light curves have 

point-to-point scatters of 0.50 mmag for the first transit and 0.59 mmag for the second. 
J> I We construct detailed signal to noise calculations for defocussed photometry, and apply 

• '"j ■ them to our observations. We model the light curves with the jktebop code and 

' combine the results with tabulated predictions from theoretical stellar evolutionary 

models to derive the physical properties of the WASP-5 system. We find that the 
planet has a mass of M b = 1-637 ± 0.075 ± 0.033 M Jup , a radius of R h = 1.171 ± 
0.056 ± 0.012 Rjup, a large surface gravity of 5b = 29.6 ± 2.8ms -2 and a density of 
Ph = 1-02 ± 0.14 ± 0.01 /9j U p (statistical and systematic uncertainties). The planet's 
high equilibrium temperature of T oq = 1732 ± 80 K, makes it a good candidate for 
detecting secondary eclipses. 

Key words: stars: planetary systems — stars: individual: WASP-5 — stars: bi- 
naries: eclipsing — methods: data analysis — methods: observational — techniques: 
photometric 
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1 INTRODUCTION 

The d iscovery of the first extrasolar planet l|Mavor fc Quelo3 
Il995h opened a new field of astronomical research. We now 
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know of over 300 extrasolar planet J], which encourages sta - 
tistical studies of these objects (e.g. lUdrv fc Santos! 120071 ). 
However, most of these discoveries have been made via mea- 
surements of the radial velocities of their parent stars, which 
yield relatively little information about individual planets 
(only the orbital parameters and lower limits on the planet's 
mass). 

The discovery of transiting extrasolar planets (here- 
after TEPs) in t he year 1999 l|Charbonneau et al.l [20001 ; 
iHenrv et alJl200d ) has provided a new window on the prop- 
erties of extrasolar gas-giant planets, making it possible to 
derive the physical properties of planets and their parent 
stars (with a little help from stellar theory). These analy- 
ses require high-precision photometry of planetary transits 
in order to measure the radii of the components and the 
inclination of their orbital plane with respect to Earth. 

Good light curves have been obtained for a substan- 
tial number of the roughly fifty known TEPs, but the pro- 
cess remains fraught with difficulties. Ground-based studies 
are subject to complications arising due to scintillation, at- 
mospheric effects which change throughout transit observa- 
tions, telescope tracking errors, and problems with satura- 
tion and flat-fielding of the ubiquitous charge-coupled de- 
vice (CCD) detectors. Space-based facilities, however, can 
be subject to major cost and calibration issues. 

We have therefore decided to experiment with heavy 
defocussing of ground-based telescopes in order to nullify as 
many as possible of the above effects. The idea is to choose 
relatively long integration times (several minutes) and dis- 
perse the resulting large numbers of photons over many 
CCD pixels. This has the huge advantage that flat-fielding 
errors can be averaged down by orders of magnitude com- 
pared to focussed observations, and also means that normal 
changes in atmospheric seeing are irrelevant. Nonlinear pixel 
responses should be corrected by calibrating the response of 
the CCD or avoided by sticking to count levels where such 
effects are low. The longer integration times and large point 
spread functions (PSFs) mean that the sky background level 
is much higher than for standard approaches, but signal-to- 
noise calculations (Section[6]l show that this is unimportant 
in many cases. 

In this work we present high-precision photometry of 
WASP-5, a V = 12.3 star which w as discovered to har- 
bour a TEP by I Anderson et aL I (|2008l ). We used the Danish 
1.54 m telescope at ESO La Silla, Chile, as part of the 2008 
MiNDSTEp campaigrfl and the DFOSC imaging camera. 
It was suspected that DFOSC would not be suited to these 
observations, as focal-reducing instruments can suffer inter- 
nal reflections, but this does not appear to be a problem. 
We observed two transit events and achieved observational 
scatters of only 0.50 mmag and 0.59 mmag per point in the 
final light curves. In the following sections we discuss our ob- 
servations and data reduction procedure s, and then analyse 
the d ata using the methods detailed by ISouthworthl (2008, 
2009). Our results are fully homogeneous with those of the 



fourteen TEPs analysed in these papers, and represent the 
highest-precision measurements of the physical properties of 
the star and planet in the WASP-5 system. 

Whilst it is unusual to break the 1 mmag barrier in 
astronomical differe ntial photometry, our observations do 
not set any records. iGilliland et al.l (|l993T ) obtained obser- 
vational scatters of 0.25 mmag per minute in time-series ob- 
servations of stars in the open cluster M67. These authors 
used a set of 4m-class telescopes, allowing them to achieve 
a lower Poisson and scintillation noise than ourselves, and 
CCD imagers with much finer pixel scale, which meant much 
less defocussing and so a lower sky background. 

A number of researchers have used some telescope defo- 
cussing to improve the photometric accuracy of light curves 
of TEPs, in many cases broadening the PSFs by only a few 
pixel s. A partial list includ e s studies of GJ 436 ( Gillon et al.l 



20071; iDemorv et alj 120071 ; lAlonso et al" 



1 See http://exoplanet.eu/ for a list of known extrasolar plan- 
ets. 

2 Information on the MiNDSTEp (Microlensing Network for the 
Detection of Small Terrestrial Exoplanets) campaign can be found 
at www.mindstep-science.org/ 



-. 12008). H P 17156 
(iGillon et al. 2008a), HP 189733 i|Winn et al.ll2007ah , HAT- 
P-l (|Winn et al.l l2007b|). Hubble Space Telescope observa- 
tions of HP 189 733~JSwain et ail 12008) ■ and the NASA 
EPOXI mission (Christiansen ct al. 2008J). The discovery 
paper of WASP-5 jAnderson et al.l 120081 ) presented a light 
curve from a 1.2 m telescope which achieved a scatter of 
about 1 mmag usi ng this techn i que. Su bstantial defocussing 
was also used by IGillon et all ([2008b) in the course of ob- 
taining a light curve of WASP-5 using the VLT, but the 
resulting data suffer systematic effects which were not re- 
movable. These were attributed to deformation of the pri- 
mary mirror, as the active optics had to be turned off 
to permit defocussing. Finally, PSF broadening has been 
performed bot h by moving the telescope during individ- 
ual exposures l|Bakos et alj|2003 ) and by shuffli ng charge at 
will round an orthogonal- transfer CCP device (|Tonrv et al.l 
120051 ; I Johnson et al.[|2008T l. 



2 OBSERVATIONS 

We observed two transits of WASP-5 using the 1.54 m Pan- 
islJ3 Telescope at ESO La Silla and the PFOSC focal- 
reducing imager. This setup yielded a full field of view of 
13.7' x 13.7' and a plate scale of 0.39" pixel -1 . For the second 
set of observations the CCP was windowed down to approx- 
imately 2000x800 pixels to decrease the readout time from 
91 s to 39 s. All observations were done through the Cousins 
R filter, which allowed a higher count rate than Johnson V 
for red stars but did not suffer the fringing effects visible 
with PFOSC and the Cousins I filter. An observing log is 
given in TableQ] The observing conditions and sky trans- 
parency were good throughout both nights. 

The exposure time was set at 120 s, to allow a large 
amount of light to be collected in each observation whilst 
avoiding undersampling the light variations through tran- 
sit (see Section|6] for detailed signal to noise calculations) . 
We then adjusted the amount of defocussing until the peak 
counts per pixel from WASP-5 were roughly 25 000 above 
the sky background. This resulted in a doughnut-shaped 



3 Information on the 1.54 m Danish Tele- 
scope and DFOSC can be found at 
http : //www. eso . org/sci/f acilities/lasilla/telescopes/dlp5/ 
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Table 1. Log of the observations presented in this work. 



Date 


Start time 


End time 


Number of 


Exposure 


Airmass 




CUT) 


(UT) 


observations 


time (s) 




2008 08 29 


05:37 


10:08 


74 


120.0 


1.03-1.54 , 


2008 09 21 


01:06 


05:32 


101 


120.0 


1.03-1.40 




Figure 1. Surface plot of the PSF of WASP-5 in an image taken 
at random from the observing sequence on the night of 2008 Au- 
gust 29th. The x and y axes are in pixels. The lowest and highest 
counts are 623 and 28429 ADUs, respectively, and the z axis is 
on a linear scale. 



PSF with a diameter of about 40 pixels (16") and contain- 
ing roughly 10 7 counts from WASP-5. Once the amount of 
defocussing was settled on, it was not changed until the end 
of the observing sequence. An example PSF of WASP-5 is 
shown in Fig.[T] 

A few images were also taken with the telescope prop- 
erly focussed, and were used to verify that there were no 
faint stars within the PSF of WASP-5 which might dilute 
the transit depth. We found that the closest detectable star 
to WASP-5 is at a distance of 44" (113 pixels), so the edge 
of its PSF is separated from that of WASP-5 by over 70 
pixels in our defocussed images. The closest stars of similar 
to or greater brightness than WASP-5 lie over 4' away. We 
conclude that no stars interfere with the PSF of WASP-5. 



3 DATA REDUCTION 

Our data reduction pipeline is written in the IDI0 program- 
ming languag e and uses the daophot photometry package 
(|Stetsorll987l ) (distributed as part of the astrolie0 library) 
to perform aperture photometry. The telescope was autogu- 
ided thoughout our observations so it was possible to fix the 
positions of the circular apertures for each dataset. We ex- 
perimented with different aperture radii, finding that this 
had a small effect on the scatter in the photometry but no 



4 The acronym IDL stands for Interactive Data Language and 
is a trademark of ITT Visual Information Solutions. For further 
details see http://www.ittvis.com/ProductServices/IDL.aspx 

5 The ASTROLIB subroutine library is distributed by NASA. For 
further details see http://idlastro.gsfc.nasa.gov/ 



detectable effect on the shape of the resulting transits. We 
chose the apertures which yielded the most precise photom- 
etry, with radii of 28 pixels for the PSF and 38 and 60 pixels 
for the inner and outer edges of the sky regions, respectively. 
Thus a PSF aperture covered 2463 pixels and a sky region 
covered 6773 pixels. 

We used these aperture sizes and the astrolib/aper 
routine to obtain aperture photometry of all stars in our 
images with sufficient counts to be useful. The number of 
comparison stars was either eight or nine for these obser- 
vations. The differential-magnitude light curves of all these 
stars were checked for short-timescale variability, and none 
was found. Most of the comparison star light curves have 
slow variations in their apparent brightness, attributable to 
atmospheric effects, so variability on timescales greater than 
about 5 hr would not be detectable. 

To remove these slow variations in the differential- 
magnitude light curves we implemented an algorithm to op- 
timise the selection of comparison stars. We defined time 
intervals outside transit and used the idl/ amoeba routine 
to minimise the square of the magnitude values within these 
intervals, resulting in light curves which were normalised to 
zero differential magnitude. The parameters of the minimi- 
sation were N coefficients of a polynomial versus time, and 
weights of each comparison star relative to the first com- 
parison star. All good comparison stars were combined into 
one ensemble by weighted flux summation. For each dataset 
we used N = 2, which removes overall variations without 
deforming the shape of the transit. We found that it was 
necessary to have N > 1 because we have too few compari- 
son stars to create an ensemble with a low Poisson noise and 
the same photometric colour as the target star. We found 
that the best-fitting comparison star weights were usually 
close to unity, as expected for observations dominated by 
Poisson noise. 

Finally, we experimented with the inclusion of bias and 
flat-field corrections to the science images. Master bias and 
flat field images were created by median-combining sets of 
bias images and dome flats. We found that subtraction of a 
master bias image had a negligible effect on the photometry, 
as the bias level of the DFOSC CCD is well-behaved and ex- 
hibits minimal changes which are simply subsumed into the 
estimation of the sky background. Including a flat-fielding 
correction, however, gives a small but noticeable decrease in 
the scatter of the resulting light curves. The flat-fielding is 
relatively unimportant in this case because the stellar PSFs 
were kept on the same pixels by autoguiding the telescope. 
We have applied bias and flat-field corrections to generate 
our final light curves, which are shown in Fig. [2] and repro- 
duced in Tabled 



4 LIGHT CURVE MODELLING 



ISouthworthl |200g| . 12009) presented an homogeneous analy- 
sis of fourteen well-observed TEPs. We have used the same 
methods here, so our findings are fully compatible with the 
results in these papers. The analysis splits naturally into 
two steps, the firs t of which is modelling the light curves 
l|Southworthll200Sh . and the second of which is including ad- 
ditional observed quantities t o determine the p hysical prop- 
erties of the star and planet (|Southworthll2009l ') . 
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Figure 2. Final light curves of WASP-5 from our two nights of 
observations. The observational error bars were too conservative 
so have been shrunk to match the scatter of the data. 



Table 2. Excerpts of the light curve of WASP-5. The full dataset 
will be made available in the electronic version of this paper and 
at the CDS. 



HJD 


Differential magnitude 


Uncertainty 


2454707.739980 


-0.000722 


0.000512 


2454707.742816 


-0.000580 


0.000513 


2454707.746034 


0.000004 


0.000512 


2454707.748360 


0.000616 


0.000512 


2454707.750814 


0.000870 


0.000513 


2454730.554371 


-0.000038 


0.000630 


2454730.556153 


0.000084 


0.000629 


2454730.557982 


-0.000347 


0.000628 


2454730.559811 


-0.000330 


0.000627 


2454730.561651 


-0.000274 


0.000625 



The transit light c urves of WASP-5 were m odelled us- 
ing the JKTEBOllfl code |Southworth et a l. 2004a. bj), which is 
based on the ebop progr am originally developed for eclips- 
ing binary star syste ms (jPopper fc Etzell Il98ll : lEtzell Il98ll : 
iNelson fe DavTI 19721 ) . The parameters of the fit include the 
fractional radii of the star and planet, which are directly de- 
pendent on the shape of the light curve and are defined to 



6 JKTEBOP is written in FORTRAN77 and the source code is avail- 
able at http://www.astro.keele.ac.uk/~jkt/ 



be ta = — and rb = — where Ra and Rh are the absolute 
radii of the components and a is the orbital semimajor axis. 
The fractional radii are actually parameterised as ta + rb 
and k = ^ in jktebop, as these are only weakly corre- 
lated with each other (|Southwortbs 2008). The other 'shape' 
parameters are the orbital inclination, i, and the limb dark- 
ening coefficients of the star, ua and va (see below). The 
mass ratio of the system is also a parameter as it affects the 
shapes of the components. We adopted a value of 0.00155, 
and verified that large changes in this number had a negli- 
gible effect on our results. 

The main analysis of our photometry was undertaken 
before any transit timings were available in the literature. 
We therefore allowed the midpoints of the two eclipses to 
float freely, by including the orbital period and eclipse mid- 
point as fitted parameters, to guard against possible period. 
We have also modelled the two light curves individually 
to measure two times of minimum light: 2454707.82388 ± 
0.00029 and 2454730.62115 ± 0.00032 (HJD). During this 
process we found that the error bars in our data (which 
come from the idl/aper routine) are slightly conservative 
and result in reduced \ 2 values of Xred =0.95 and 0.75 for 
the two datasets. We therefore multiplied the error bars by 
yj Xrcd De f ore combinin g the data from th e two nights into 
one light curve (see e.g. lBruntt et al]|2006l ). 

Our two times of minimum light have been au g mented 
with the four transit times given by iGillon et~al] (l2008bh 
and used to determine a linear ephemeris. lGillon et al.l found 
that a straight line was not a good representation of the data 
(Xred = 5.7). With our additional data we can associate 
this discrepancy with the transit time from their FTS light 
curve, which is earlier by 3<r that predicted by our linear 
ephemeris. Further data are needed to decide whether the 
cause is astrophysical or observational. Using all times of 
minimum light, we find 

T = HJD2 454 375.62494(24) + 1.6284246(13) x E 

where quantities in parentheses denote the uncertainty in 
the fin al digit of th e prece ding number. 

In lSouthworthl (2008) it was found that the treatment 
of the stellar limb darkening (LD) requires careful thought. 
LD is included in jktebop as a choice of several paramet- 
ric 'laws' for which coefficients must be specified. The linear 
law is inadequate for high-precision observations, whereas 
more sophisticated nonlinear laws suffer from strong corre- 
lations between t heir coefficients (So uthworth et aT] |2007a; 
ISouthworthll2008T ) . Furthermore, theoretically calculated LD 
coefficients depend on stellar model atmospheres and do not 
match the highest-quality observations. The solution is to 
adopt nonlinear laws but to fix one of the coefficients at a 
reasonable value. An alternative would be to re-parameterise 
the two-coefficient laws so the coefficients have weaker cor- 
relations, but this has little effect on the other parameters 
of the fit. 

We have incorporated LD using the linear, quadratic, 
square-root, logarithmic or cubic laws. We present one set of 
solutions where both the linear (ma) and nonlinear (va) LD 
coefficients are fixed to theoretical values (Table|3]). We also 
give a set of solutions where ua is a fitted parameter and va 
is fixed to theoretical values (Table|3|). In the second case, va 
is perturbed by ±0.1 on a flat distribution during the error 
analysis. The theoretical LD coefficients are taken to be the 
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Table 3. Parameters of the JKTEBOP best fits of the light curve of WASP-5 for different LD laws with the coefficients fixed at theoretically 
predicted values. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. 
Results were not calculated using the cubic LD law because theoretical LD coefficients are not available. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


»"A + »"b 


0.2036 ±0.0056 


0.2009 ±0.0057 


0.2085 ±0.0050 


0.2024 ±0.0057 


k 


0.1103 ±0.0010 


0.1096 ±0.0010 


0.1119 ±0.0007 


0.1100 ±0.0009 


i (deg.) 


86.17±0.84 


86.66 ±0.96 


85.37 ±0.60 


86.37 ±0.90 


"A 


0.60 fixed 


0.40 fixed 


0.20 fixed 


0.70 fixed 






0.30 fixed 


0.50 fixed 


0.20 fixed 


VA 


0.1834 ±0.0048 


0.1811 ±0.0051 


0.1875 ±0.0044 


0.1824 ±0.0050 


rb 


0.02022 ±0.00071 


0.01984 ±0.00070 


0.02098 ±0.00061 


0.02006 ±0.00071 


a (mmag) 


0.5696 


0.5593 


0.5586 


0.5602 



Table 4. Parameters of the JKTEBOP best fits of the light curve of WASP-5 for different LD laws with the linear coefficients included 
as fitted parameters and the nonlinear coefficients fixed at theoretically predicted values which are perturbed in the Monte Carlo error 
analysis. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


+ n> 


0.2097 ±0.0052 


0.2034 ±0.0062 


0.2058 ±0.0053 


0.2054 ±0.0056 


0.2059 ±0.0060 


k 


0.1125 ±0.0009 


0.1102 ±0.0013 


0.1112±0.0011 


0.1110±0.0011 


0.1114 ±0.0012 


i (deg.) 


85.16 ±0.64 


86.22 ±1.02 


85.77 ±0.75 


85.83 ±0.77 


85.73 ±0.81 


u A 


0.508 ±0.031 


0.378 ±0.048 


0.227 ±0.049 


0.659 ±0.058 


0.495 ±0.037 


VA 




0.30 perturbed 


0.50 perturbed 


0.20 perturbed 


0.10 perturbed 


VA 


0.1885 ±0.0046 


0.1832 ±0.0055 


0.1852 ±0.0047 


0.1849 ±0.0049 


0.1852 ±0.0052 


n> 


0.02119 ±0.00067 


0.02019 ±0.00083 


0.02058 ± 0.00070 


0.02053 ±0.00071 


0.02063 ±0.00077 


cr (mmag) 


0.5576 


0.5579 


0.5578 


0.5579 


0.5578 



tainties to the paramete rs found from the l ight curve, fo l- 
lowing the procedures of ISouthworth et all (|2004bl . l2005bl) . 
Each fit was followed by 1000 Monte Carlo simulations 
and la error bars were calculated from the scatter in the 
simulation results (Tables [3] and [4}. This does not ac- 
count for any systematic errors which might be in our light 
curves, so this possibilit y was investigated u sing the residual- 
permutation approach l|jenkins et al.|[2~002l ) and found to be 
marginally significant. The residuals in Fig.0 suggest there 
are minor systematic effects, in agreement with the residual- 
permutation results. 



From Tables [3] and |4] it can be seen that the solutions 
with LD coefficients fixed at theoretical values have slightly 
larger residuals than those where the linear LD coefficient 
is a parameter of the fit. We have therefore used the lat- 
ter solution s for our final resu lts, combined following the 
methods of ISouthworth! |2008l ) and incorporating the sys- 
tematic uncertainties found with the residual-bead method. 
The final values of the photometric parameters are given 
in TableO They are in agreement with and more precise 
th an the values found in the t wo previous studies of WASP- 
5, lAnderson et~aT] (|2008h and ICillon et al] l|2008bl ), which 
used the same photome tric data. Includ i ng the stellar veloc- 
ity amplitude given bv lAnderson et al] (|2008h . we measure 
the surface gravity of the pla net to be <?b = 29.6 ± 2.8 ms -2 
(see ISouthworth et al.ll2007bl ). 




-0.06 



-0.04 



-0.02 0.00 0.02 
Orbital phase 



0.04 



0.06 



Figure 3. Phased light curve of the two transits of WASP-5 
compared to the best fit found using JKTEBOP and the quadratic 
LD law with both LD coefficients included as fitted parameters. 
The residuals of the fit are plotted at the bottom of the figure, 
offset from zero. 



approximate mean of val ues found by bilinea r interpolation 
within the tabulations oflVan Hammd (Il993l l . IClaret] (|2000l. 
l2004al ) and lClaret fc Hauschildt] (|2003l ). 

We have used Monte Carlo simulations to assign uncer- 
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Table 5. Final parameters of the fit to the light curve of W ASP-5, based on the indiv idual solutions given in Table[J] The orbital 
ephemeris was determined separately (Section|4j. The results of lAnderson et aT, ] (20081) and lGillon et al.l ll2008rj) . which are based on the 
same photometric data, are included for comparison, and are without error bars if they were not directly quoted quantities. To is the 
reference time of mid-transit. 



This work 



Anderson et al. f2008 



Gillon et al. 



Orbital period (d) 
T (HJD) 

ta + rb 
k 

i (°) 



1.6284246 ± 0.0000013 
2454375.62494 ± 0.00024 
0.2052 ±0.0068 
0.1110 ±0.0014 

85.8 ±1.1 
0.1847 ±0.0061 
0.02050 ±0.00091 



1.6284296 1 



h0. 0000048 
-0.0000037 

2454375. 62466+ ' 00026 



-0.00025 



0.197 



1092+ - 0030 
u.iuyz_ Q 0017 

85.0 to 90.0 

0.178 

0.0194 



1.6284279 
2454373.99598 



+0.0000022 
0.000004!) 
+0.00025 
-0.00019 

0.198 



n insfi+ - 0010 
0.1086_ 0013 

3+2.8 
-0.7 



86. 9: 
0.179 
0.0194 



5 PHYSICAL PROPERTIES OF WASP-5 

The photometric parameters determined above have been 
used to infer the physical proper ties of the WASP- 5 sys- 
tem using the method outlined in ISouthworthl |2009). It is 
not possible to calculate the properties directly from obser- 
vations, as one piece of information is missing and must 
be filled in by additional constraints. This approach re- 
quires as input the measured values of ta, rb, i, orbital 
period, and the stellar effective temperature, metallicity 
and velocity amplitude. For the latter three quantities we 
adop ted T eB = 5 7 00 ± 1 00 K and [f] = 0.09 ± 0.09 dex 
from iGillon et all d2008bh and K A = 278 ± 8kms _1 from 
I Anderson et al.h 20081 ) . We then interpolated within the tab- 
ulated predictions of several sets of theoretical stellar mod- 
els to provide the optimum fit t o the input quantities. The 
stella r mo dels used were Claret JClaretJl2004bl. 120051. 12006] 
2007]) . Y 2 (|Demaraue et al.ll2004h and Padova jGirardi et all 



2000). The uncertainties in the input parameters were prop- 



agated into th e output physical properties using a perturba - 
tion analysis l|Southworthl l20od ; ISouthworth et all l2005al ) . 
which allows a detailed error budget to be obtained for each 
output quantity. 

The use of theoretical stellar model predictions means 
that the calculated physical properties of the WASP-5 sys- 
tem are subject to systematic errors caused by any de- 
viations of the models from reality. Unfortunately, stel- 
lar models are persistently unable to match the directly 
measured physical properties of low-mass eclipsing binary 
star syst ems, although this effect is l ikely due to magnetic 
activity (lHoxidll973l ; IClausenl 1 19981 ; iTorres fc RibasH2002l ; 
iRibas et al j2008j T A lower limit can be set on the systematic 
errors by considering the variation in results from different 
sets of stellar models. This is only a lower limit because inde- 
pendent models have some physical ingredients in common, 
for example opacity tables. A probable upper limit can be 
set by avoiding stellar models and instead constraining the 
properties using an empirical main-seq uence mass-radius re- 
lation for the star (|Southworthll2009l ). Tabled] contains the 
physical properties of WASP-5 calculated using the three 
different sets of stellar models and with the mass-radius re- 
lation. 

The properties calculated using stellar models are in 
good agreement with each other, whereas those determined 
using the mass-radius relation are offset to larger values. 
The difference is most significant for the orbital semimajor 



Table 7. Detailed error budget for the calculation of the physical 
properties of the WASP-5 system from the light curve parameters, 
stellar velocity amplitude, and the predictions of the Claret stellar 
models. Each number in the table is the fractional contribution 
to the final uncertainty of an output parameter from the error 
bar of an input parameter. The final uncertainty for each output 
parameter (not given) is the quadrature sum of the individual 
contributions from each input parameter. P represents the orbital 
period. 



Output 




Input parameter 






parameter 


P K A 


i 


rA 


rb 




m 


Age 






0.341 




0.753 


0.556 


a 


0.001 




0.075 




0.760 


0.644 


M A 






0.075 




0.760 


0.645 


Ra 






0.871 




0.374 


0.317 


l°g3A 






0.967 




0.195 


0.166 


PA 






1.000 








M b 


0.627 


0.031 


0.059 




0.592 


0.502 


R h 






0.028 


0.927 


0.284 


0.241 


9b 


0.307 


0.015 




0.951 






Pb 


0.208 


0.010 


0.010 


0.969 


0.099 


0.083 



axis (a) and the stellar mass (Ma) as these quantities are 
quite affected by systematics but almost unaffected by ob- 
servational uncertainties. The error budget for the output 
parameters is given in Table[7] and shows that more precise 
stellar T g and [^] measurements would be useful, as well 
as further high-precision photometry. 

Table[8] gives the final physical properties for WASP-5, 
with statistical errors obtained from the perturbation anal- 
ysis and systematic errors from consideration of the intera- 
greement between results calculated using the three different 
sets of stellar evolutionary models. We have adopted the val- 
ues of the properties obtained using the Claret models, in or- 
der to retain homogeneity with the analyse s of ISouthworth 
(2009). A comparison w i th the results of lAnderson et al. 
(2008) and IGillon et all (|2008bl l shows that all quantities 
agree within the quoted uncertainties. The main difference 
is that our more extensive photometry yields a larger radius 
and thus a lower density for the planet; this arises from the 
light curve analysis via the parameter rb- 
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Table 6. Physical properties for WASP-5, derived using either an empirical stellar mass-radius relation or the predictions of different sets 
of stellar evolutionary models. The stellar mass, radius, surface gravity and density are denoted by Ma, Ra, l°g<?A an d PAi respectively. 







Mass-radius 


Padova models 


Y 2 models 


Claret models 


a 


(AU) 


0.02777 ±0.00078 


0.02702 ±0.00046 


0.02728 ± 0.00043 


0.02729 ±0.00049 


M A 
R A 
l°g9A 


(M Q ) 

(Rq) 
[cgs] 


1.076 ±0.091 
1.103 ±0.060 
4.385 ±0.024 


0.991 ±0.051 
1.073 ±0.039 
4.373 ±0.030 


1.020 ±0.048 
1.083 ±0.040 
4.377 ±0.029 


1.021 ±0.055 
1.084 ±0.039 
4.377 ±0.030 


PA 


(P©) 


0.803 ±0.080 


0.803 ±0.080 


0.803 ±0.080 


0.803 ±0.080 


M b 

Rb 
9b 
Pb 
Age 


(M Jup ) 
(Rjup) 

(ms-2) 

(PJup) 

(Gyr) 


1.70±0.11 
1.191 ±0.063 
29.6±2.7 
1.00 ±0.14 


1.604 ±0.073 
1.159 ±0.055 
29.6 ±2.8 
1.03 ±0.14 


1.635 ±0.070 
1.170 ±0.055 
29.6 ±2.8 
1.02 ±0.14 

5 9+ 2 ' 2 


1.637 ±0.075 
1.171 ±0.056 
29.6 ±2.8 
1.02±0.14 

6.8± 3 ;° 



Table 8. Fi nal physical propertie s for WASP-5. The first e rror bars are statistical and the second are systematic. The corresponding 
results from lAnderson et alJ l|2008h and iGillon et al ] d2008bh have been included for comparison. 



Final result (this work) Anderson et al. (2008) Gillon et al 



a 

M A 
Ra 

logSA 
PA 

M b 
^b 
9b 
Pb 



(AU) 

(M ) 

(Rq) 

[cgs] 

(pq) 

(M Jup ) 

( Rjup) 

( ms -2 ) 

(PJup) 



0.02729 ± 0.00049 ± 0.00027 



0.02683 



1.021 ± 
1.084± 
4.377 ± 
0.803 ± 
1.637 ± 
1.171 ± 
29.6 ± 
1.02 ± 



0.055 

0.040 

0.030 

0.080 

0.075 

0.056 

2.8 

0.14 



±0.030 
±0.011 
±0.004 
±0.000 
±0.033 
±0.012 

±0.01 



+0.00088 
0.00075 
+0.099 



972 

u ' 3 ' -0.079 

1 026+ - 073 
i.uzd_ 044 

4 4QQ+0.039 

0.90 

1 58+° 13 

1 090+ 094 

30.5+ 32 
i 99 +0.i9 

1 - zz -0.24 



0267+ 0012 

U - UZD( -0.0008 

96+ 13 
UMD -0.09 

-, qoq+0.056 

4 395+ - 043 
* ,o3J -0 .040 

0.88 ±0.12 

1 co+0.13 

30.5+ 4 ° 

O+0.26 



1.23 



-0.16 



6 WHAT IS THE OPTIMAL EXPOSURE 
TIME? 

Having presented and interpreted our observations, we now 
present detailed signal to noise (S/N) calculations which can 
be used to plan similar projects. The output quantities in- 
clude the amount of defocussing required, and the total noise 
to signal ratio in millimagnitudes. The tractability of these 
calculations demands several approximations, so the results 
should not be overinterpreted. The idl code used for the 
calculations can be obtained from the first author. 

The important point here is that the quantity of interest 
for a given observing sequence is the achievable S/N per unit 
time rather than per observation. This is limited in particu- 
lar by CCD readout for short exposures and sky background 
noise for long integrations. There will be an optimum value 
between these two extremes, which depends on the target 
and observing facilities. 

We start with count rates for the target (overall) and for 
the sky background (per pixel) gathered from the SlGNAlQ 
code. These are appropriate for the Isaac Newton Telescope 
and Wide Field Camera, so must be scaled for the differ- 
ence in telescope aperture (1.5 m versus 2.5 m) and plate 
scale (0.39" pixel -1 instead of 0.33" pixel - ). We adopt the 
default values for sky brightness in dark, grey and bright 
time; this approximation can be improved as and when ad- 



7 Information on the Isaac Newton Group's SIGNAL code can be 
found at http://catserver.ing.iac.es/signal/ 



ditional information is available. The count rates are speci- 
fied in electrons, so the CCD gain level is taken into account 
with the input parameters rather than as part of the S/N 
calculations. 

For a given t cxp (in seconds) we calculate the total elec- 
trons obtained from the target, and the associated Poisson 
noise, to be 



Starget - 
A'target 



^expCta 



(1) 

(2) 



where t oxp is the exposure time and Ctargot is the electrons 
per second from the target. In order to estimate the amount 
of defocussing required (the actual shape of the PSF does not 
enter these calculations), we distribute these equally over 
n p ix pixels by specifying a maximum number of electrons 
per pixel from the target and the background, m to tai. Fig.[T] 
shows that this approach is imperfect so we also specify the 
maximum electrons per pixel from the target alone, m t arget- 
This allows us to limit the range of counts within a struc- 
tured PSF, and thus avoid saturating the most illuminated 
pixels. 

If the sky background is low then we need only to ac- 
count for the limit set by m t argct in working out how many 
pixels to put the target electrons into: 

/ \ ^target . . 

n p ix(m t argot ) = (3) 

^Titarget 

If the sky background is high then we must increase the 
number of pixels containing counts from the target to avoid 
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broaching the m to tai limit. In this case the number of pixels 
in the PSF is: 



Ptotal — Ototal 



n P ix(mtotai) 



"St a 



m-total ~ Csky^cxp 



(4) 



where C a k y is the electrons per second per pixel from the 
sky background. For a specific situation, n p i x is taken to be 
the larger of the two alternatives. 

Now we have found n p i x we can calculate the total sky 
level in the PSF, and the ensuing Poisson noise, in one ob- 
servation: 



5sky — £cxp7ipixC^sky 



-^sky ^/^cxp^pixCsky 

The contribution from CCD detector readout noise is: 

-^Vr-on — ^ron \/"mx 



(5) 
(6) 

(7) 



where n ron is the readout noise per pixel in electrons. 

The treatment of flat-field noise is more difficult, as the 
PSFs of one star in different observations will have many 
pixels in common whose flat-field noise will therefore cancel 
out. We approximate this noise contribution using /flat, the 
flat-field noise per pixel (expressed as a fraction of the elec- 
trons in a pixel) , and apply it to a number rifl at pixels where 
"-flat ^ ripix- The flat-field noise is then 



, i ^target . \ n flat 

/flat I — 1" Osky 



(8) 



We can now express these noise levels as fractions of the 
total electrons from the star in magnitude units: 



7 sky 



""flat 



= -2.5108! 
-2.51og 10 
-2.5 log 10 
-2.51og 10 



Sta 



Stars 



sky 



^target 
Starget Ar ( 



Sta 



•Sta 



A^fla 



Sta 



(9) 
(10) 

(11) 
(12) 



The nois e due to scintillation c an be found using the r elation 
iven by iDravins et alj fl998) and originally due to lYound 
1963): 



■^"scint 



: 0.004D" 2/3 X 7/4 e"' l/H (2t exp ) 



-1/2 



(13) 



where D is the telescope aperture (m), X is the airmass, h 
is the altitude of the telescope (m), and H = 8000 m is the 
scale height of the atmosphere. Note that scintillation varies 
seasonally, and the exponent of the airmass term actually 
varies between 1.5 and 2.0 depending on t he direction of 
the w ind relative to the telescope pointing (jPravins et al.l 
1998). These caveats illustrate the approximate nature of 
S/N calculations. 

We now have five noise contributions which can be 
added in quadrature to reveal the total noise for one mag- 
nitude measurement in one observation: 



O-total = y' of arget + of ky + (Tron + <Tfl at + ^int (14) 

The ratio of noise to signal per unit time is then 



(15) 



where d roa dout is the total dead time per observation (nor- 
mally dominated by the CCD readout time). In reality there 
will be additional noise in differential photometry due to 
measurement of the magnitude of the comparison star. We 
do not include this here as its contribution to the total noise 
can be driven to a much smaller value than for other sources 
by using many comparison stars. 

Using the above formulae we have performed S/N cal- 
culations for the observations presented in this work. We 
adopted m ta rgct = 10 000, M tota i = 40 000, n ron = 4.0, X = 
1.2, and neglected flat-fielding noise. We set d roa dout = 91 s 
to represent the first of our two observing sequences. The 
predicted Ototai for t oxp = 120 s is 0.54 mmag (dark time), as 
compared to the 0.50 mmag we actually achieved. This good 
agreement could be further improved by tinkering with the 
input parameters to the S/N calculations. 

The resulting curves of noise level per observation and 
per minute are plotted as a function of i exp in Fig. [4] The 
'sweet spots' representing the lowest achievable noise level 
per minute (indicated by filled circles) occur at £cxp — 565 s 
(dark time), 310 s (grey) and 163 s (bright). Specifying a 
shorter readout time or fainter target star would shift the 
sweet spots towards smaller t cxp . The abrupt change of slope 
visible in Fig.|4]at t cxp = 740 s (bright time) occurs when the 
sky background becomes so high that the starlight has to 
be spread over extra pixels to avoid transgressing the mtotai 
limit. Using these S/N calculations we find that, except for 
a target star which is very faint, defocussing is always the 
best approach. 

The exposure times yielding the lowest noise per minute 
for the current situation are too long to sample the light 
variations through a transit of WASP-5 adequately, which 
is why used t eX p ^ 120 s. This is partly due to the long 
readout time of the DFOSC CCD: adopting drcadout = 39 s 
(appropriate for our second transit sequence) shortens the 
optimal exposure times to 371 s (dark) , 203 s (grey) and 107 s 
(bright). 



7 DISCUSSION AND CONCLUSIONS 

We have presented high-precision light curves of two transits 
of the transiting extrasolar planetary system WASP-5, ob- 
tained by heavily defocussing the Danish 1.5 m telescope at 
ESO La Silla. The defocussing caused the telescope PSF to 
take an approximately annular shape of diameter 40 pixels 
(16"). Data reduction was undertaken using standard aper- 
ture photometry techniques. Debiassing and flat-fielding the 
raw CCD images improved the results but only by a small 
amount, probably due to the use of autoguiding to confine 
the stellar PSFs to the same pixels on each image. Our final 
light curves have a scatter around the best-fitting model of 
0.50 mmag and 0.59 mmag for the two transits, observed in 
dark and grey time, respectively. 

We have presented detailed signal to noise calculations 
for defocussed photometry, taking into account Poisson noise 
from the target and background, readout noise, scintilla- 
tion, and flat-fielding noise. Application of the resulting for- 
mulae to the current case yields a predicted noise level of 
0.54 mmag (dark time), which is in good agreement with the 
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Figure 4. Predicted noise levels for the observations presented in 
this work, but as a function of exposure time. The dotted curves 
show the predicted noise level per observation for dark, grey and 
bright time (lower to upper curves). The solid curves show the 
predicted noise level per minute of telescope time (dark, grey, 
bright), and the optimum exposure times are indicated with filled 
circles. 



actual values. We use these calculations to assess what expo- 
sure times would give the best overall signal to noise per unit 
time, finding 565s (dark time) to 163s (bright time). Such 
large values arise partly because of the long readout time of 
the DFOSC CCD, and are impractical because they would 
not properly sample the light variations through transit. 

We find that telescope defocussing is a viable way of ob- 
taining high-quality light curves of point sources over short 
periods of time. Its main advantages are that flat-fielding 
noise is averaged down and a smaller proportion of the ob- 
serving time spent reading out the CCD. The scintillation 
and Poisson noise levels per unit time are ultimately lim- 
ited by the telescope aperture rather than the observational 
strategy. One possible disadvantage of our approach is the 
use of a focal-reducing imager. These can be prone to inter- 
nal reflections, which may have slightly increased the noise 
levels we find. One definite advantage, though, is that the 
telescope was equatorially mounted, so the target and com- 
parison stars could be kept to the same light path whilst 
observing a transit. It is not possible to do this with altitude- 
azimuth telescopes, even if an image derotator is used, mak- 
ing such observations liable to problems with systematic er- 
rors. 

The resulting light curves were mod elled with th e jkte - 
bop code using the approach outlined bv lSouthworthl {2008). 
Limb darkening of the star was dealt with using five differ- 
ent LD laws, with similar results when one or both of the 
coefficients were included as parameters of the fit. We then 
adopted the observed Jeff, [H] and velocity am plitude of the 
star (|Anderson et alJl2008tlGillon et all feoOSb) and interpo- 
lated within the tabulated predictions of three different sets 
of theoretical stellar evolutionary models to derive the phys - 
ical properties of the WASP-5 system jSouthworthl [2009V 
This approach has allowed us to assign both statistical and 
systematic errors to all output quantities. Our results are 
fully consistent with the homogeneous analyses of 14 tran- 



siting extrasolar planetary systems presented bv lSouthworthl 

(|20odl20ogh . 

We find that the planet WASP-5 b has a density equal 
to that of Jupiter, and that the system has an age of ap- 
proximately 6 Gyr. The theoretical models presented by 
iFortnev et aD \200H ) can match WASP-5 b's mass and ra- 
dius without requiring a heavy-element core. WASP-5 is also 
in line with the trend of planets with sh ort orbital periods 
to hav e relatively hi gh surface gravitie s l)Southworth et aD 
l2007bl ) and masses (|Mazeh et al.ll2005l ). We calculate that 
its equilibrium temperature is T cq = 1732 ± 80 K, which may 
be high enough for WASP-5 b to fall into the 'pM class' of 
planets with large atmo spheric opacities du e to oxides of 
titaniu m and vanadium (IFortnev et al.ll2008l ). Its ISafronovl 
ljl972t ) number is 9 = 0.0755 ± 0.0096 , which makes it a 
Class II planet a ccording to the stud y of lHansen fc Barmanl 
(120071 1. However. ISouthwort hi l|2009l ) found that the division 
between Class I and Class II objects is under threat from 
the steady addition of new discoveries to the catalogue of 
known transiting extrasolar planets, and may not be sta- 
tistically significant. The high equilibrium temperature of 
WASP-5 b makes the system a good ca ndidate for detecting 
the secondary eclipse a t phase 0.5 (e.g . Deming et alJfeOOa ; 
lHarrington et al I l2007l ; iKnutson et all 120071 ) . which would 
help our understanding of the physics of the atmospheres of 
gas giants. Such observations, as well as more precise mea- 
surements of the stellar velocity amplitude, T c g and 
are encouraged. 



"Ml 
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